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In a recent article, we constructed a hierarchy Bl of outer boundary conditions for Einstein's 
field equations with the property that, for a spherical outer boundary, it is perfectly absorbing 
for linearized gravitational radiation up to a given angular momentum number L. In this article, 
we generalize £>2 so that it can be applied to fairly general foliations of spacetime by space-like 
hypersurfaces and general outer boundary shapes and further, we improve B2 in two steps: (i) we 
give a local boundary condition C2 which is perfectly absorbing including first order contributions 

^ — , in 2M / R of curvature corrections for quadrupolar waves (where M is the mass of the spacetime and 

R is a typical radius of the outer boundary) and which significantly reduces spurious reflections due 
to backscatter, and (ii) we give a non-local boundary condition T>2 which is exact when first order 

£\J , corrections in 2M/R for both curvature and backscatter are considered, for quadrupolar radiation. 

PACS numbers: 04.20.-q, 04.25.-g, 04.25.Dm 
O 1 I. INTRODUCTION AND MAIN RESULTS 

m ; 

I Most often, numerical simulations of black hole spacetimes involve replacing the physical unbounded domain with 
. a compact domain £ with artificial outer boundary <9£. In order to obtain a well-defined Cauchy evolution, boundary 
Q\ " conditions need to be specified at 9E. Ideally, these conditions should be formulated in such a way that the artificial 
boundary is completely transparent to the physical problem on the unbounded domain. However, in reality, one aims 
for boundary conditions which (i) form a well-posed initial boundary value problem (IBVP) and (ii) insure that very 
little spurious reflection of gravitational radiation occurs from the outer boundary. Boundaries which meet these two 
criteria are called absorbing. 

The formulation of absorbing outer boundary conditions in general relativity has always been a challenging problem, 
primarily because constraint violations can enter the numerical domain via the outer boundary, the definition of in- 
and outgoing radiation is ambiguous, and the evolution of the geometry of the outer boundary is not known a priori. 
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Strides towards the development of absorbing outer boundaries which both preserve the constraints and control the 
^Jy gravitational radiation by freezing the Newman-Penrose Weyl scalar ^0 to its initial value have been made in Refs. 
Ill) 0, B BL H @- For other approaches whichpreserve the constraints exactly but control the radiation at the 
boundary through different means, see Refs. 0, II, S 03 • In [ll|, boundary conditions which are perfectly absorbing 
for quadrupolar solutions of the flat wave equation are analyzed and numerically implemented via spectral methods, 
and proposed to be used in a constrained evolution scheme of Einstein's field equatio ns Other methods avoid 

introducing an artificial outer boundary altogether by compactifyin g sp atial infinity 13, [l4| or by making use of 
hyperboloidal slices and compactifying null infinity (see, for instance |l5l. flrj fl7l|). 

Recently, a hierarchy Bl of local boundary conditions for \T/ has been proposed (see Eq. (128) of [TH) which, for 
L > 2, improves the freezing ^0 boundary condition given in [l], 0, S 0j S Q m that it is perfectly absorbing on a 
spherical boundary for linearized gravitational perturbations with angular momentum number £ < L. If the outer 
boundary is placed in the wave zone (meaning kR 3> 1, where R is a characteristic radius of the outer boundary and 
k is the characteristic wave number) and the solution is smooth enough, the few lower £ multipoles will dominate. 
Therefore, the boundary condition Bl with L = 2, 3 or 4 should yield very little spurious reflections in linearized 
gravity. When considering the full nonlinear problem, if the outer boundary is not only in the wave zone but also far 
away from the strong field region (meaning R^> M where M is the total mass of the system, that is, the Arnowitt- 
Deser-Misner (ADM) mass), one expects the spacetime to be described by a linear perturbation about Minkowski 
spacetime. Therefore, in the nonlinear case, the boundary conditions Bl for L — 2, 3 or 4 should also yield small 
spurious reflections provided that kR 3> 1 and R^> M. 

In this article, we formulate Bl to be applicable to more general slicings than in [1 81 ] ; furthermore, we construct 
Bl in such a way that the outer boundary is not restricted to be an approximate metric sphere. The approach we 
adopt is to describe a general four dimensional metric near the outer boundary of the computational domain as the 
Schwarzschild metric plus a perturbation. We perturb Schwarzschild, rather than Minkowski, spacetime in order to 
incorporate leading order effects in M/R into our boundary conditions. For an asymptotically flat spacetime, the 
M/R correction terms to Minkowski spacetime are given by the monopole contribution of the Schwarzschild metric. 
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When these terms are taken into account, a gravitational wave solution on Minkowski spacetime acquires two types 
of correction terms. The first is the curvature correction term which obeys Hugyens' principle and the second is a 
fast decaying term which violates Hugyens' principle and describes the backscatter off of curvature. 

A method which provides exact outer boundary conditions for linearized waves propagating on a Schwarzschild 
background has been recently presented in [19.l20l.l2lj. However, in view of the nonlinear theory, there is no advantage 
to obtaining a boundary condition which takes into account the exact form of the Schwarzschild metric beyond the 
order of M/R, since the full spacetime agrees with the Schwarzschild metric only up to order M/R. If second order 
effects are to be taken into account, then linear effects in J/R 2 (where J is the total angular momentum of the system) 
and quadratic effects in M/R from the full metric should also be included. 

Our boundary conditions Bl are formulated as follows. We assume that in the vicinity of the outer boundary, 
the spacetime metric is the Schwarzschild metric of mass M plus a small perturbation of it. More precisely, we 
assume that the spacetime manifold near the outer boundary can be written locally as a product of a two-dimensional 
manifold M and a unit two-sphere S 2 such that if (x a ) denote local coordinates on M and (x A ) — ip) denote the 
usual angular coordinates on S 2 , the components of the full metric g^ ul1 ^ satisfy the following conditions: 



(full) ~ I „ , 
9a b ~ 9ab\ < 1, 



\ 9 { Ib U) \ 



«1, \g { lB l) /r 2 -9AB\«l, 



(1) 



where g = dfl 2 + sin $d<p 2 is the standard metric on S 2 , r is the areal radius, and g is a pseudo-Riemannian metric 
on M representing the radial part of the background Schwarzschild metric. In standard Schwarzschild coordinates 
(t, r), we have g = —N(r)dt 2 + N(r)~ 1 dr 2 , where N = 1 — 2M/r. However, we do not assume that (x a ) are 
Schwarzschild coordinates: with respect to the coordinates (x a ) on M, the components of g are given by the coordinate 
transformation of the (t, r)-components of g which minimizes \g„t u — g a b\- If the boundary is spherical, a natural 
representation of the spacetime manifold in the neighborhood of the outer boundary as M x S 2 exists, and our hope is 
that one can enforce clever boundary conditions on the gauge degrees of freedom guaranteeing that conditions (TT]) are 

satisfied. In this case, a method for computing r and g from the full metric is the following (2^]: if a a b = 9^^ an d 
Iab = 9ab denote the induced metrics on M and S* 2 , respectively, then r is defined by Anr 2 = J s2 y / py|sin(^)(i^(i</j 
and g :— J s2 a si^^dddcp. If the boundary is not spherical, g and r need to be computed by other methods: a 
possible approach is to use a minimization technique in order to optimize the inequalities in ([T|). 

Having identified M, g and r, one introduces the coordinates t and r* on M. which correspond to the Schwarzschild 
time and tortoise coordinates, respectively, but which are defined in a completely geometric way (see Section [TTJ) so 
that different observers on M agree on their definition. Consequently, the boundary condition Bl does not depend 
on the foliation of M . It reads 



Br 



dt[r 2 (d t + d r ,)] 



L-l 



5 *o) 



= 0, 



(2) 



where the Weyl scalar "Jo is defined in terms of the Newman-Penrose null tetrad given in Sect. HH Notice that Bi is 
the freezing boundary condition. 

After formulating Bl, we concentrate on improving B2 by reducing spurious reflections off the outer boundary due 
to curvature corrections and backscatter. We accomplish this by computing in- and outgoing quadrupolar solutions 
for even and odd-parity perturbations of the Schwarzschild metric g + r 2 g, to first order in 2M/R. Using these 
approximate solutions, we compute reflection coefficients which quantify the amount of spurious reflections. When 
calculating the reflection coefficients, we assume for simplicity that the outer boundary is an approximate metric 
sphere with constant area 4irR 2 . This allows us to compare the amount of reflections resulting from our improved 
boundary conditions with those resulting from B\ , for a spherical outer boundary. Our boundary condition Z? 2 already 
gives a 1.2{M / R)(kR)~ 1 times reduction in the amount of spurious reflections compared to B\ for kR > 1.05. 1 We 
find we can improve B2, however, by applying the operator in Eq. ([2]) to r 5 iV _1 ^'o rather than to r 5 ^ (recall 
N = 1 - 2M /r) giving 



C-2 



9t(5 t + 9 r J(r 5 AT- 1 * )| as = 



(3) 



For a spherical outer boundary, C2 reduces spurious reflections due to curvature corrections and backscatter by a factor 
of (15M/2R)(kR)~ 2 compared to Bi, for kR > 1. We have not calculated the reflection coefficients for boundary 
shapes other than spherical; however, we stress that our boundary conditions Bl and C2 can be applied to general 



1 In our previous article [Til ], we found only an (M/R) times reduction in the amount of spurious reflections compared to B\ for kR > 1.05. 
This discrepancy is due to the fact that the calculations in 1 1 SH were performed with a slightly different version B2 of 62- 
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outer boundary shapes. Furthermore, we show that in general, C2 is perfectly absorbing to first order in 2M/R if 
backscatter is neglected but the curvature correction terms are taken into account. 

Finally, we improve C2 so that it is perfectly absorbing up to order 2M / R including backscatter for quadrupolar 
gravitational radiation. This boundary condition reads 



r(d t + d r J(r 5 N- l *o) + 



30M 



t/2r 



r 5 ^ (t-2ry,r,-d,tp)dy 

(i + y) 7 



G, 



(4) 



where G is some appropriate boundary data which can be computed from the portion of the initial data which is 
exterior to the outer boundary (see Eq. (|38p below). In particular, G — if the initial data vanish in the region 
exterior the the outer boundary. 

Our article is organized as follows. In Sect. [Til first we consider an arbitrary spherically symmetric background 
manifold of the form Ai x S 2 and show how, in vacuum, the coordinates t and r* can be defined in a geometric 
way. This definition is important as it makes it possible to apply our boundary conditions to quite general foliations 
of the spacetime and to spherical as well as non-spherical outer boundary shapes. Next, using the gauge-invariant 
perturbation formalism presented in [23], we derive the linearized Weyl scalars S^o and 8^4 when there are both even 
and odd-parity perturbations on a Schwarzschild background in arbitrary spherically symmetric coordinates. These 
scalars are invariant with respect to infinitesimal coordinate transformations and tetrad rotations, and we use S^q to 
construct our new boundary conditions. In Sect. Hill we compute the in- and outgoing quadrupolar solutions to the 
Regge- Wheeler and Zerilli equations up to order 2M/R. Here, we find the interesting result that for the perturbed 
Weyl scalar 5^o, the ingoing solutions are identical for both the even and odd-parity sectors, whereas the outgoing 
solutions for these sectors differ by a sign in a first order correction term. In Sect. IIV1 we discuss the boundary 
conditions B2 and C2 and calculate the corresponding reflection coefficients up to first order corrections in 2M/R. For 
B 2 , the reflection coefficient decays as (2M/R)(kR)~ 5 and for C 2 , it decays as (2M/R)(kR)~ 6 (for kR > 1). Finally, 
we construct the boundary condition T>2, which is perfectly absorbing including backscatter and curvature correction 
terms, to first order in 2M/R. Conclusions are drawn in Sect. fVl where we also present the boundary condition Cl 
[2^ | , a generalization of C2 derived from the results in [25j . In the appendix, we analyze the stability of the boundary 
condition 2?2 for a related model problem. 



II. LINEAR PERTURBATIONS OF A SPHERICALLY SYMMETRIC VACUUM SPACETIME 

In this section, we first consider an arbitrary spherically symmetric background spacetime which can be represented 
as the product of a two-manifold M. and a two-sphere S 2 and show how the coordinates t and r* on M can be 
defined geometrically, provided that the background spacetime satisfies Einstein's vacuum equations. We then use 
the perturbation formalism in [23|, which is a generalization of the Regge- Wheeler [26| and Zerilli [27j equations 
to arbitrary spherically symmetric coordinates, to derive formulas for calculating the linearized Weyl scalars S^/q, 
<5\E f 4, and Im(5 1 i>2)- As constructed in this section, these Weyl scalars are gauge-invariants in the sense that they are 
invariant with respect to infinitesimal coordinate transformations and tetrad rotations. 

The generalized Regge- Wheeler and Zerilli equations found in [23| have the form of wave equations with potentials 

for the gauge-invariant scalars and 3>«i which describe, respectively, odd and even-parity linear fluctuations 

of the background geometry with angular momentum numbers Im. In [23j, it was also shown that the complete set 

of metric perturbations in the Regge- Wheeler gauge can be reconstructed from the scalars Q^J 1 and $11, without 
solving additional differential equations. In particular, the linearized Weyl scalars 5'^q 1 ^4 an d \xa{§^%), which are 
gauge-invariant quantities, are entirely determined by these scalars. 



A. Background equations 

We consider a spherically symmetric background of the 2 + 2 form M x S 2 with metric 

g = g a bdx a dx b + r 2 c/ab dx A dx B , (5) 

where g = di) 2 + sin 2 §d(p 2 is the standard metric on S 2 and g and r denote the metric tensor and a positive function, 
respectively, defined on the two-dimensional pseudo-Riemannian orbit space M. In what follows, lower-case Latin 
indices refer to coordinates on M, while capital Latin indices refer to the coordinates $ and <p on S 2 . The ah 
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components of the Einstein tensor belonging to the metric ([5]) are 

Gab = --V V 6 r + ^g ab ferAr + N - l) , (6) 

where V denotes the covariant derivative compatible with g, A = V a V a and N — g(dr, dr) = g ab V a r ■ V r. A 
coordinate-invariant definition of the ADM mass is given by 

M= r -{l-N). (7) 



In vacuum, Eq. ^ yields 



M 

V a V b r = -^~g ab . (8) 



Eqs. (O and ([8|) imply that M is constant. 

For the following, we assume that N > 0. Since in this article we are interested in the asymptotic regime of the 
Schwarzschild spacetime, where N is close to one, this assumption poses no problems. The two vector fields 

t a = _£^v 6 r, r a = V°r (9) 

on M, where e a b denotes the natural volume element on M, are orthogonal to each other and are linearly independent 
since r a r a = —t a t a = N > 0. As a consequence of Eq. ([5]), \7 a tb — Ms ao /r 2 ; hence, t a is a Killing vector field. 
Furthermore, it follows from ([8]) that t a and r a commute. Therefore, it is possible to introduce coordinates (i, r*) on 
M such that t a d a = dt and r a d a — d rt . The coordinates t and r* correspond to the standard Schwarzschild time and 
tortoise coordinates, respectively. For later use, we also introduce the two null vector fields 

l a d a = -^(t a + r a )d a = -j={d t + d r „), (10) 

k a d a = -l=(t a -r a )d a =-^=(d t -d r J, (11) 

satisfying g a b — —lakb — k a l b and e a b — lak b — k a h- They can be completed to form a Newman-Penrose null 
tetrad {l a , k a ,m B ,fh B }, where m B is a complex null vector field and fh B its complex conjugate such that r 2 gAB = 
uiArriB + friA'mB and the volume element i belonging to g is given by r 2 EAB = i(mArriB — m^ms). The Weyl scalars 
are defined as 

*o = C a0li J a m^m u , 

*2 = C a p ia/ l a mPmf l lt v , 
*3 = C af3 ^l a k"k^rh 1 ', 

where Cap^ denotes the Weyl tensor. For a spherically symmetric spacetime of the form ([5]) with the tetrad choice 
(jlOllip . all the Weyl scalars vanish except for 

M 

*2 = T • 



B. Perturbation equations 

Next, let us consider a generic linear perturbation of the background metric: 

5 g a b = L ab , 8g a B = QaB , SgAB = r 2 K A B , 

where the quantities L a b, QaB and Kab depend on the coordinates (x a ,x B ). With respect to an infinitesimal 
coordinate transformation Sx a t— > Sx a + £ a , 5x B i— > Sx B + i} B , these quantities transform according to 

Lab H-» Lab + 2V (a £ o) , (12) 
QaB QaB + VsCa + r 2 V a 77B , (13) 

2 

K AB K AB + - §ab r a £ a + 2V (a t? b) , (14) 
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where V denotes the covariant derivative compatible with g and where t)a = gABV B - A lengthy but straightforward 
calculation yields the following expressions for the perturbed Weyl scalars S"i>o and S^4'- 



S*o = V abCD ri b m c m D 



= V abCD k a k b fh c m D , 



where V abC D is given by 



VabCD = g (-V c V d L q6 + V c V (a Q b)D + V D V {a Q b)c ~ V {a r 2 V b) K CD 



if 



(15) 



(16) 



and where the superscript tf refers to the trace- free part with respect to the metric gcD- It can be verified readily that 
VabCD is invariant with respect to the infinitesimal coordinate transformations given in Eqs. (|12I13I14|) . Furthermore, 
since S^o and 6^4 do not depend on the perturbed tetrad components, they are also invariant with respect to 
infinitesimal tetrad rotations. In fact, this invariance of S^o and 5^ 4 holds for any type D background metric [281 ]. 

Since the background metric is spherically symmetric, it is convenient to expand the fields in spherical harmonics. 
In the resulting perturbation equations, pieces belonging to different angular momentum numbers I and m decouple. 
Thus, it is sufficient to consider one fixed value of £ and m at a time. The decomposition of the fields L ab , Q a B and 
Kab into tensor spherical harmonics reads 



L ab — H ab Y, 
QaB = Q a V B Y + h a S B , 



(17) 



K AB = Kg AB Y + G 



VaV b 



tf 



Y + 2kV {A S B 



where Y — Y ip) denotes the standard spherical harmonics and S A = £a V bY . With respect to this, we obtain 



V a 



bCD 



-H 



[inv) 
ah 



tf 



Y + W {a h^W (c S 



(C D D) , 



(18) 



where H^™^ and h^ 71 ^ are the gauge-invariant amplitudes defined in [23[ which, in the Regge- Wheeler gauge, reduce 
to the quantities H ab and h b , respectively. These gauge- invariants were first given in [29j . For vacuum perturbations, 
these gauge- invariant quantities are entirely determined by two scalar fields <E>( + - ) and $( _ ) according to 



H 



(inv) 



A + a 



VaV fc -ig ab A) r(A + a )$ (+) 



Ainv) _ ~ ™ 
"ft — £afcV 



r$(~) 



(19) 

(20) 



where A = (£— i)(t+2) and a = 6M/r. The two scalar fields <&( +) and <&( _ ) satisfy the Zerilli and Regge-Wheeler 
[2^ | equations, respectively: 



where the potentials are given by 

V+(r) = 
V-(r) = 



-V Q V a + V±(r) $ (±) =0, 

A 2 r 2 [£(£ + l)r + 6M] + 36M 2 (Ar + 2M) 

(Ar + 6A/) 2 r 3 
£{£ + 1) 6M 



(21) 



r" r" 

Taking everything together, we obtain the following expressions for the perturbed Weyl scalars S^q and 8^4: 



1 



1 



r z I A + ao 



S^A = ~77 



n fo v a v t 



r 2 I A + a 



A a fc b V a V fc 



r(A + a )$ (+) 
r(A + a )$ (+) 



il a l°V a V b 
f ik a k b V a V 



m c m D V c V D Y, 



m m VcVdi, 



where we have defined 



and where the scalar fields $( ± ) satisfy Eq. (|2T 



Finally, we compute the imaginary part of the perturbed Weyl scalar 6^2- One first finds 

Im(<S* 2 ) = -lm(W abC Dl a k b m c m D ) , 



(22) 
(23) 

(24) 
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where the quantity 



WabCD — 



V C V [Q 



Qb]i 



VbVu 



Qb]c 



(25) 



is invariant with respect to the infinitesimal coordinate transformations given in Eqs. (|12I13I14|) . Next, application 
of the decomposition (fTT|) yields 



WabCD = r 2 V [a 



h 



(inv) 
b] 



Using Eq. fl20j! and 2r 2 m c fh D V \ c Sd] = »9 VcV D F, we then obtain 



Im(5*2) = - 



1(1+1) 
4r 



V°V a $ (_) - — [ 1 - — J $ ( ~ } 



Y. 



Finally, we substitute the Regge- Wheeler equation ([2Tj) into the above expression to get 



Im(J* 2 ) 



{£-!)£{£ +!){£ + 2) 
4r 3 



$(-)y. 



(26) 



Therefore, the gauge-invariant scalar $(~~) obeying the Regg e- Wheeler equation can be interpreted as the radial part 
of the imaginary part of the perturbed Weyl scalar 8^2 [la. |30||. Comparison of Eq. ([26)) with the quantity tj> = r 2 ho 
defined in [l8| gives the relation 2(f> — {I — !)£{£ + 1){£ + 2)&~\ It would be nice to know what the geometric 
interpretation of the gauge-invariant scalar obeying the Zerilli equation is. 



III. IN- AND OUTGOING QUADRUPOLAR WAVE SOLUTIONS 

In this section, we compute in- and outgoing solutions to the Zerilli and Regge- Wheeler equations. Since we 
are interested only in the asymptotic regime r ~ R, with R large compared to M, we compute these solutions to 
first order in 2M/R. A method for performing this task is described in [l8|]. We start by computing the outgoing 
solutions. For this purpose, we introduce the outgoing Eddington-Finkelstein coordinates (r,p), which are related to 
the geometrically defined coordinate t and the areal radius r of Sect. Ill Al through 



r = t + r — r* 



p = r, 



where the tortoise coordinate r* is defined by r*(r) = r — 4M+2Mlog(r/2M— 1). In these coordinates, the two-metric 
reads 

g = -dr 2 + dp 2 + - — idr - dp) 2 . 
P 

With respect to this metric, the Zerilli and Regge- Wheeler equations (f2Tj) can be written as 



dl - 3 2 



£(£+!) 



$(±) 



2M 
P 



(d T + d p ) 2 --(d T +d p )~^ + of 2}l 
P P 



V R 



(27) 



where a+ = (A + 4)/ A and a_ = 1. In order to compute the solution to first order in 2M/R, we make the ansatz 
(omitting the superscript (±) on <!>) 

^. , _ . . 2M . 
$(r, p) = $ (r, p) + -g-fll(r, p), 

where $0 is a solution of the Regge- Wheeler-Zerilli equations (f2Tj) for 2M/R = 0, and where the function g\ satisfies 
the equation 



9i(t,p) = - — 
P L 



1 



3a± 



(dr+d p ) 2 --(d T +d p )-^- 

p p 



$0 



(28) 
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In the following, we compute $ for quadrupolar radiation (1 = 2). In this case, the outgoing flat background solution 
$0 is given by 

<&o(t, p) = U<V(p - t) - luM (p - t) + \u{p - r), 

P P 

where U is a smooth function on the real axis with fc'th derivative . In the following, we also assume that U is 
bounded and that it vanishes for negative enough values of its argument. The corresponding solution for g± is given 
by 



gi(r, P) - 2 f )R U{p - r) + ^fr (1) (P - r) + | / K-A r. ,, .nVU)dr. 



p-T 



where the integral kernel K% reads 



(29) 



K 2 (t, p,x) = -- 1 h~ 4 + 2w- 3 + 2w- 2 ] 



" 2p 



X > p — T, 



and satisfies 



p 1 
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K 2 (t,p,x) = 0, K 2 (T,p,p- t) = — 



The solution <?i is not unique; one can add an arbitrary homogeneous solution of Eq. ([28)1 to it. However, g\ is 
uniquely characterized by the following conditions: 



lim gi(const. + p, p) = (<?i vanishes at future null infinity), 

p — >oo 

lim gi(const. — p 1 p) — (<?i vanishes at past null infinity), 

p — >oo 



lim gi(r, const.) 



(g± vanishes at future time- like infinity). 



Notice that the third requirement is necessary in order to exclude solutions of the form gi(r, p) = (ciT + c )p~ 2 (where 
Co and c\ are some constants) of the zero mass Regge- Wheeler- Zerilli equations (f2Tj) , which vanish at both future and 
past null infinity. 

Summarizing, we have the following outgoing approximate solutions of the Regge- Wheeler- Zerilli equations for 
1 = 2: 



<p { ±\t,r) = u^{n-t)--U^'{n-t) + ^u ± {n-t) 



r(2) 



2M 



%M , (r.-*) + *£^C«'.-.) + ! 



K 2 (r — r* + 1, r, x)U±(x)dx 



r»-t 



O 



where a + = 2, a_ = 1. Since the Regge- Wheeler-Zerilli equations (|2ip are time-symmetric, corresponding ingoing 
solutions are obtained from this by merely flipping the sign of t: 



$^(t,r) = vi 2) (r*+t)--^i 1) (r*+t) + - 3 2^±(^+*) 



2Af 



^Kf) (r . +() + 3£^>5v ± (r. +t ) + f 



if 2 (r — r* — i, r, x)V±(a;)dx 



/2M 
V 1T 



Using equation (j2"2")) and the identity 



which holds for the Schwarzschild background, we obtain the following expressions for the perturbed Weyl scalar S^>q: 



6V Q 



1 



(30) 
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where 



3 1, s 2M 



-2l7±(r* -t)T^£ ) (r*-t) + ^ / A-( I +/y)f"i (/■ - / -2/-;/),/// 



/2MV 



\ R J 



(31) 



r 4 iV 

2M 
r 



l/ ± (r* + t) - 2rV^ 1} (r* + t) + 2r^ ZJ (r* + 1) - ^^(r* + 1) + ^rVf(r. + *) 



„3t/(3). 



„4t/(4), 



1 /" V"±(r* +t + 2ry)dy 



yf 



o 



(2M 



(32) 



with k(w) = 5w~ 6 + 4u>~ 5 + 3w^ i + 2w~ 3 + ; 



In the odd-parity case, these expressions agree with the expressions 



obtained in [18| except that {/_ and V_ are rescaled by a factor of 12. This factor is explained by Eq. (f2"fj)) . Notice 
that the expressions in the even (+) and odd (— ) case are exactly identical except for the sign in the third term of 
V>o/' Finally, we notice that d^o/* decays as r~ 5 along the outgoing future-directed null geodesies while 8^>o\ decays 
as r _1 along the outgoing past-directed null geodesies, as expected from the peeling theorem 31]. 

The first two terms inside the square brackets in the expressions for ipQ^} and describe the curvature corrections, 
while the integral terms describe the backscatter. Notice that, by construction, the solution given in Eq. (|3"2"|) does 
not yield any radiation at future null infinity and hence the integral in that expression is acausal in the sense that it 
involves integrating over the infinite future of V±. However, since we are only needing the zeroth order contribution 
in M / R of 3/>o\ in this article, this acausality does not introduce any problems here. 



IV. ABSORBING BOUNDARY CONDITIONS 



In this section, we focus on the boundary condition £>2 which, as formulated in this article, is applicable to general 
outer boundary shapes and fairly general foliations of the spacetime. We improve Bi to give two new boundary 
conditions, C2 and T>2, which reduce the amount of spurious reflections due to curvature and backscatter off of 
curvature. C2 is obtained by modifying B2 so that it is perfectly absorbing for quadrupolar gravitational radiation up 
to order 2M/R curvature (but not backscatter) terms. We find that the reflection coefficient for C2 is smaller than 
that for B2 by a factor of the order of kR, where k is the characteristic wave number of the gravitational wave. Finally, 
we construct T>2, which is perfectly absorbing up to order 2M/R in both curvature and backscatter for quadrupolar 
waves. By nature of the backscatter effect, T> 2 is non-local: as we will see, it involves an integral over the past 
boundary. 



A. An estimate for the amount of spurious reflections 



Here, we estimate the amount of spurious reflections for quadrupolar waves for the boundary condition B2, and 
compare it to the corresponding result for the freezing ^0 boundary condition B\ . In order to calculate the reflection 
coefficients, we make the following three assumptions: 

(i) the boundary is an approximate metric two-sphere with constant area 47ri? 2 , 

(ii) the outer boundary lies in the weak field regime, i.e. R 3> 2M, 

(iii) if k is the characteristic wavenumber of the wave, Mk <C 1 . 

With assumptions (i) and (ii), we can linearize the field equations about the exterior Schwarzschild spacetime of mass 
M near the outer boundary, where M is the total mass of the system, and assume that the outer boundary is located 
at constant areal radius r = R. In- and outgoing quadrupolar wave solutions near the outer boundary are described 
by the expressions found in the previous section. The general solution to the IBVP with the boundary condition B2 
consists of a superposition of an in- and a outgoing solution: 



V-'o = ipoy + 7 ipo\ , 



(33) 
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with ipoy an d V'ov gi yen by Eqs. (f3Tj) and (|32|) . respectively, and where 7 is an amplitude reflection coefficient. To 
determine 7, we assume monochromatic radiation, 



ik(r,+t) 



(34) 

with k a given wave number satisfying < k <C 1/M, and plug the ansatz (f3"3")l into the boundary condition i3 2 - First, 
we find (the result in the even- and odd-parity case is exactly identical, so we omit the superscript (±)) 



r(9t + $..)(r 4 W0(*,r) 



6Af 
r 



C/(n - 1) - 15 



f7(r* - t + 2ry)dy 



(i + y) 7 



r(d t + d r ,)(r^)(t,r) = ir 5 V&(r* + 1) + O {^j 
Using the monochromatic ansatz (|34[) and imposing #2, we then find 



O 



/2M 



7 = 7 kR, 



2M 
~R~ 



with the integrals 



C n (z) 



3M e 2lkR * 
' 2iR (kR) 5 

cos(2zy) 



[1 - l5C 7 (kR) - 15iS 7 (kR)} + O 



flM 



(i + yy 



dy, 



S n {z) 



sin(2zy) 

(i + y) n 



dy, 



n > 2, 



which have the asymptotic expansions (2z) 2 C n (z) = n[l — (n + l)(n + 2)/(2z) 2 + 0(z 4 )], 2zS n (z) = 1 — n(n + 
l)/(2z) 2 + 0(z~ A ) for z > 1. Finally, we obtain 



7 fcK, ■ 



2Af 

IT 



2Af 



E{kR) + O 



/2M 



(35) 



where 



£(z) = A [(1 - 15C 7 (z)) 2 + (155 7 (z)) 2 



1/2 



Because of assumption (iii), we may neglect the quadratic corrections in 2M/R. (If kM is of the order of unity or 
larger, powers of kR which are multiplied by (2M/R) 2 might be comparable in size to or larger than terms of the 
form M/R times unity. For example, see the exact outgoing solutions on a Schwarzschild background obtained in 
[25|].) For kR ^> 1, the reflection coefficient goes as (2M / 'R)(kR)~ 5 . Comparing this result with the corresponding 
result for B\ , for which the reflection coefficient uM 



\l2{kR)\ = 



l--(kR) 6 + -(kRf 



-1/2 



O 



(f) 



decays as (fci?) -4 for kR 3> 1, we see that we gain a factor of (2M / R)(kR)~ 1 in the decay rate. The function 
(kR)E(kR)/\j2(kR)\ is plotted in Figure [T] This plot, together with the asymptotic expansion 2zE(z)/\j2(z)\ — 
1 + 7/(8z 2 ) + 0(z~ 4 ), suggests that for kR > 1.05, this function does not exceed 0.6. Therefore, if kR > 1.05, the 
boundary condition B2 gives a reflection coefficient which is smaller by a factor of 1.2M / (kR 2 ) than the one for By. 

One can do better: applying the operator r(dt + d r ,) to the function r 5 N^ ll i' Q instead of r 5 ^/ gives the boundary 
condition C 2 given in Eq. ([3]). Since in this case, outgoing quadrupolar wave solutions satisfy 



r(d t +d r ,)(r 4 N~ 1 ^)(t,r) 



90M f U(r* - t + 2ry)dy 



o 



(2M 
\R 



(36) 



the term 6MU(r* — t)/r is eliminated. Physically, this means that the boundary condition C2 takes care of the 
curvature correction terms and is perfectly absorbing if the backscatter off the curvature is neglected. The function 
E(z) in Eq. (|35p is now replaced by the function 



F(z) 



45 

4^5 



[C 7 (zf + S 7 (z) 



21 1 /2 
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so the new reflection coefficient decays as (2M/R)(kR) 6 as kR — > oo. Therefore, we gain a further power of (kR) 1 
in the decay rate! The function (kR) 2 F{kR)/\^ 2 (kR)\ is plotted in Figure[TJ Together with the asymptotic expansion 
Az 2 F(z)/\ l2 (z)\ = 15[1 - 71/(8z 2 ) + 0(z~ 4 )], this plot suggests that for kR > 1, {kR) 2 F{kR)/\ l2 {kR)\ does not 
exceed 15/4. Therefore, if kR > 1, the improved boundary condition C 2 gives a reflection coefficient which is smaller 
by a factor of 15M/(2k 2 R 3 ) than the one for B x . 
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FIG. 1: Comparison of the functions E(kR), E(kR), and F(kR) (which are the reflection coefficients times R/(2M) for the 
boundary conditions B2, B2, and C2) divided by the reflection coefficient |72(fc-R) for the freezing- boundary condition Si, 
versus kR. Notice that the reflection coefficients are scaled by factors of kR for the sake of the presentation. 

For comparison purposes, we plot the function E(z)/\'y2{z)\ in Figure[U where the function E(z) was computed in 
[H[ and represents the reflection coefficient resulting from the boundary condition B 2 ■ d t (d t + d r )(r 5l $o) | r=7? = 
(the same as B 2 but with d rt replaced by d r ). We have checked that - as with the functions E(z) and F(z) - the 
function E(z) is exactly the same in the odd- and the even-parity sector. 



B. New boundary conditions which reduce spurious reflections due to the backscattering 

Here, we show how to obtain even better boundary conditions: ones which are perfectly absorbing up to order 2M/R 
including the backscatter for quadrupolar gravitational waves. In other words, we now construct boundary conditions 
which have the property of giving a reflection coefficient that is of the order of (2M/R) 2 for quadrupolar radiation. 

First, by comparing Eq. flM]) with ip /(t, r) = 3r~ 4 [U(r* -t) + 0(2M/R)}, we notice that the boundary condition 

00 

m + ».x^'».)U + 2£ / ^"^y^ = (37) 



is satisfied up to first order in 2M/R for the outgoing solution, the integral operator taking care of the leading order 
corrections in the backscatter. However, a problem with the boundary condition in Eq. (|37p is that it requires 
knowledge of ^0 on the whole past portion of the boundary surface. So for a Cauchy formulation which starts an 
evolution from initial data at the time slice t = 0, say, this boundary condition is impractical. In order to fix this 
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problem, we split the integral in the expression (I36[) into two parts: 



t/2R oo 

U(R*~t + 2Ry)dy _ [ U(R* - (t - 2Ry))dy f U(R* + 2Ry - t)dy 



(i + y) 7 J (i + y) 7 J (i + y) 7 

U t/2R 

For a fixed event (t, R) with t > 0, R ^> 2M, the first integral on the right-hand side can be interpreted as an integral 
over the past line (r, R), < r < t, whereas the second integral can be interpreted as an integral over the line (0, r), 
R < r < oo. With this in mind, we replace the boundary condition (|37[) with the condition: 



r(fi t + $..)(r B JV- 1 tf )| 



t/2R oo 

30M f R 5 ^ (t-2Ry,R,ti,<p)dy 30M /" s 5 * (0, s, tf, </?) 



dy, 

s={R t +2R,y-t). 



R J (i + y) 1 R J (i + y) 7 

t/2R 

(38) 

where (.). : (— oo, oo) — > (2M, oo) denotes the inverse of the transformation (.)* : (2M, oo) — * (—00,00) which maps 
r to r,. In particular, (r,)* = r and (r*). = r. Now the integral on the left-hand side of Eq. (|38[) only involves the 
past portion (t,R), < r < t, of the boundary point (t,R). The past portion (r, i?), r < 0, which is not available 
in a Cauchy evolution starting at i = 0, is replaced by an integral over the exterior r > R of the initial data. If the 
initial data are compactly supported in the interval (0, R), this integral is zero and can be discarded. If the initial 
data do not vanish for r > R, one uses the exterior part of the initial data in order to define the integral on the 
right-hand side of Eq. (|38p. which acts as a source term in the new boundary condition. The generalization of Eq. 
(|38|) to non-spherical boundaries yields the boundary condition T> 2 given in Eq. ([4]). 

Summarizing, we have constructed a new boundary condition which is perfectly absorbing up to order 2M/R in 
the backscatter for quadrupolar gravitational radiation. In the appendix, we show for a related model problem that 
this boundary condition is stable in the sense that it admits an energy estimate. 



V. CONCLUSIONS 



In this article, we have analyzed the problem of specifying absorbing outer boundary conditions in general relativity. 
To this end, we have constructed approximate solutions to the Einstein equations linearized about a Schwarzschild 
metric of mass M (where M describes the ADM mass of the total spacetime), using the generalized Regge-Wheeler- 
Zerilli formalism of Ref. [23[ . These solutions represent in- and outgoing gravitational radiation in the asymptotic 
regime 2M/r -C 1, where r denotes the areal radius. Since far enough from the strong field region, any station- 
ary asymptotically flat spacetime looks like a Schwarzschild spacetime of mass M, these solutions describe, up to 
order 2M/r, in- and outgoing gravitational radiation propagating on the asymptotic region of any such spacetime 
background. 

There is an interesting difference between the propagation of even and odd-parity outgoing gravitational linearized 
radiation on a Schwarzschild background. This difference can be obtained by computing the scalars and 
on M. in the decomposition 



l a l b V abCD = r^ +) 
k a k b V abCD = r^| +) 



tf 



tf 



Y-r^ ] V (C S D) 
Y + ryj { - ) V (c S D) 



of the gauge-invariant quantity V a bCD defined in Eq. (|16p (recall that (+) refers to even and (— ) to odd parity). Using 
Eqs. (J3TJ) , ([23]), and the identity 2r 3 fc a fc b V a Vf,(r$) = N[r 2 N- 1 {d t - <9 r J] 2 $, we obtain the following expressions for 
the outgoing solution: 



r±N 



s 2M 

U±(u-t) + 



(39) 



U±{r* - 1) - 2rf/i 1) (r st - t) + 2r z U^>{u - 1) - ~r- 3 £4° J (r* - *) + ^ 4 ^± V* - *) 



„2rr(2) 



3rr(3), 



,4rr(4), 
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2M 
r 



i)-ir 3 E/i 3) (r* 



U±(r*-t + 2ry)dy 

(i + y) 2 







(40) 



with iV = 1 — 2M/ r and fc(iu) = 5w -6 + 4w~ 5 + 3k;~ 4 + 2w~ 3 + w~ 2 . Taking into account the different normalization 
of the null tetrad vectors l a and k a , and identifying G(u) — 1 2 \U+ (r* — t) + iU- (r* — t)] with u = t — r*, expression ([40]) 
agrees precisely with the expression given in Eq. (4.26) of Ref. [25| |. up to order M/R. However, the identification 
in Ref. [2|| of G(u) in Eq. (4.18) with G(u) in Eq. (4.26) is incorrect. Expression (|59|) agrees with the one given in 



Eq. (4.18) of Ref. |25j with the identification G(u) = 3 [£/+(>* - t) + iJ7_(r* - t)} - 3M[U ( + ] (r» - t) - iU { } ] (r, - i)]/2. 
While the expressions for are exactly the same in the odd and even parity sectors, those for the quantity V'o/' 
differ by a sign in the third term. 2 For M = 0, there is no difference: even- and odd-parity waves behave in exactly 
the same way. But as soon as the curvature of the background is taken into account, there is a symmetry breaking 
between the two parity sectors. Provided sufficient accuracy is available, one should be able to detect this difference 
in actual numerical simulations. 

Using the ingoing and outgoing solutions, and the assumption that the outer boundary is an approximate metric 
sphere of constant area, we have computed the reflection coefficient which quantifies the amount of spurious radiation 
reflected into the computational domain by our boundary condition B2 given in Eq. ([SJ • While this boundary condition 
is perfectly absorbing for linearized quadrupolar waves, there are small spurious reflections when the curvature of the 
background is taken into account. However, we have shown that the reflection coefficient due to this effect is extremely 
small; it is about a factor M /(kR 2 ) smaller than the corresponding reflection coefficient for the freezing V^o boundary 
condition which was already found to be small in [l8[ . By slightly modifying B2 , we obtain a new boundary condition 
C2 which results in even less reflections. To give a specific example, if the outer boundary is spherical with areal 
radius R > 100M and quadrupolar waves with wavelength 50M are considered, the freezing boundary condition 
yields a reflection coefficient which is smaller than 6.5 x 1CU 5 . If the boundary condition C2 is used instead, the 
reflection coefficient is smaller than 3.1 x 10~ 8 . Finally, we have found a boundary condition T>2 which takes into 
account the leading order behavior of the backscatter and is perfectly absorbing up to order 2M/R for quadrupolar 
gravitational radiation. However, this boundary condition is non-local: it involves two integral operators. The first 
operator involves an integral over the past of the boundary point, the second an integral over the exterior portion of 
the initial data. It is shown in the appendix for a related model problem that such non-local boundary conditions are 
stable in that that they admit an energy estimate. 

Although our reflection coefficient calculations assume a spherical outer boundary, the formulation of our boundary 
conditions do not: all that is needed for their construction is that spacetime in the vicinity of the outer boundary can 
be represented as the product of a two-manifold M. and a two-sphere S 2 such that the assumptions |T]) hold. 

We have constructed the boundary conditions C2 and T>2 to reduce spurious reflections for quadrupolar waves. 
However, in many physically interesting scenarios such as the binary black hole problem, it is possible that when 
implementing T>2, spurious reflections from the £ = 3 octupolar contribution to the gravitational radiation could 
be as large or larger than the £ = 2 quadrupolar backscatter correction. Bardeen [24[ has generalized C2 using the 
results in [25} to give Cl , a hierarchy of local boundary conditions which is perfectly absorbing including curvature 
corrections (but neglecting backscatter) to order 2M/R for all multipoles of gravitational radiation up to a given 
angular momentum number L. He obtains 



as 



0. 



(41) 



To determine whether C3 or T>2 would more effectively reduce spurious reflections when octupolar contributions are 
taken into account, we have calculated the reflection coefficient due to backscatter of quadrupolar waves for C3: 



7 kR 



2M\ 

~r) 



45 M 
4(fci?) 5 "R 



(6C 7 (kR) ~ 7C s (kR)) 2 + (6S 7 (kR) - 7S 8 (kR)Y 
(kR) 2 + 9 



1/2 



(42) 



Using the asymptotic expansions of C n (z) and S n (z) given in Sect. IIV A[ we find that Eq. (|4*2"j) decays as 
(2M/R)(kR)~ 7 as kR — > 00. Although not calculated, the reflection coefficient due to backscatter of octupolar 



2 This result contradicts the Teukolsky-Starobinsky relation given in Eq. (41) of chapter 9 in [32l |. which does not distinguish between odd- 
and even-parity perturbations. Based on the generalized Regge-Wheeler-Zerilli formalism used here, we have re-derived the Teukolsky- 
Starobinsky identities for Schwarzschild perturbations and verified that our result is correct. In particular, we have calculated the 
constant C in Eq. (41) of chapter 9 in |32| and found that it is different in the odd- and even parity sectors. 
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waves for C3 is expected to decay at least as fast or faster. For the boundary condition T>2, we find that when oc- 
tupolar radiation is considered, the reflection coefficient is (to leading order in 2M/R) the same as that given in Eq. 
(96) of [H[ and falls off as (kR)~ 6 for large kR. We conclude that indeed C3 would be more effective than P 2 unless 
the octupolar contributions of the gravitational radiation are significantly suppressed. Using the results in [251 ] . the 
backscatter corrections to first order in 2M/R can in principle be calculated to give U3 and even T>l- 

When applied to full nonlinear formulations of Einstein's field equations, the boundary conditions Cl and T>2 on ^>o 
should be used along with constraint-preserving boundary conditions and boundary conditions on the gauge degrees of 
freedom so that the resulting initial-boundary value problem is well posed. A possible strategy for specifying boundary 
conditions on the gauge degrees of freedom is to insure that the outer boundary is a metric sphere throughout the 
evolution. This should yield a natural split of the spacetime manifold into a two-dimensional manifold M and the 
two-sphere S 2 . If the outer boundary is not a metric sphere, our boundary conditions can still be applied; however, the 
above split of the manifold first needs to be identified. Well posed initial-boundary value formulations incorporating 
the boundary conditions Cl and T>2, and their numerical implementation, will be studied in future work. 
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APPENDIX A: ON THE STABILITY OF NON-LOCAL BOUNDARY CONDITIONS 

In this appendix, we analyze the following initial-boundary value problem on the quarter space Q R :— (0, 00) x 
(-00, R), R > 0: 

u{t,x) - u"(t,x) + V(t,x)u(t,x) = F(t,x), (t,x)€to R , (Al) 

OO 

u(t, R) + u'{t, R) + jpf k R {t, s)u(s, R)ds = G(t), t > (A2) 


u(0,x) = f(x), u(0,x)=g(x), x < R. (A3) 

Here, a dot and a prime denote differentiation with respect to t and x, respectively, F is a given smooth source 
function on fl R , G is a given smooth source function on (0, 00) and / and g are smooth initial data on (— 00, i?). 
Furthermore, we assume that the potential V in Eq. (|A1|) and the integral kernel k R (., .) : (0, oo) 2 — > R in Eq. (|A2[) 
satisfy the following three conditions: 

(i) V : £Ir — ► R is measurable and bounded in the following sense: There exists a smooth and bounded function 
W : (—00, R) — ► R and a strictly positive constant p > such that p/R 2 < W{x) for all x < R, V(t, x) < W(x) 
for all (t, x) € Q,r and 

R 2 

D := sup — [W(x) - V(t, x)] < 00. 

(If V is smooth, bounded, does not depend on t and satisfies V(x) > p/R 2 for all x < R we can take W(x) = V(x) 
and D = 0.) 

(ii) kft satisfies the causality condition knit, s) = if s > t, 

(iii) There is a constant C > which is independent of R such that 

00 

J \k R (t, s)\dt < CR for all s > 0, 



00 

J \k R (t, s)\ds < CR for all t > 0. 
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By a suitable re-definition of u and F, if necessary, we can assume that the boundary condition is homogeneous, that 
is, G = 0. In this case, and under the assumptions (i),(ii) and (iii) above, we now show that a sufficiently smooth 
solution of the IBVP (|Al|A2IA3|) which vanishes for x sufficiently negative fulfills the energy estimate 



£{t) < e"*/ fl 



t R 

£{Q) + cirJ J \F{s,r)\ 2 drds 

-oo 



(A4) 



where the energy is defined by 



£{t) = - I (\ii(t,x)\ 2 + \u'(t,x)\ 2 + W{x)\u(t,x)\ 2 )dx, 



and where the dimensionless constants a and /i are given by 

1 



26 



t i = 6 + D 



C 2 



6 being a strictly positive but otherwise arbitrary constant. The estimate (|A4[) proves uniqueness and continuous 
dependence of the solution on the data. We do not prove existence of solutions here. 

In order to prove the estimate (|A4[) we first differentiate £ (t) with respect to t and use the evolution equation ()A1[) 
and assumption (i) and obtain 



dt 



£{t) 



{u(t, x)u"{t, x) + u'(t, x)u'(t, x) + [W(x) - V(t, x)} u(t, x)u(t, x) + u(t, x)F(t, x)} dx 



< u(t,R)u'(t,R) + ^^£(t) + ^ J \F(t,x)\ 2 dx, 

— OO 

where 5 > and where we have used Schwarz's inequality in the last step. Integrating over t and using the boundary 
condition f|A2|) . we obtain 

t t oo 

£(t)-£(0) < - J \u'{s,R)\ 2 ds- J J k R (s,T)u'(s,R)u(T,R)dTds 



D + 5 
R 



£(s)ds 



R 
26 



o o 

t R 



\F(s,x)\ 2 dxds. 



(A5) 



In order to estimate the second term on the right-hand side, we use assumption (iii) and Schwarz's inequality again, 
and obtain 



kR{s, t)u'(s, R)u(t, R)drds 



o o 



< / \u'(s,R)\ / Ikn^T^lkR^r^lu^R^drjds 



1/2 



1/2 



< / \u'(s,R)\ / \k R (s,T)\dr 



\k R {s,T)\\u(T,R)\ 2 dT\ ds 



1/2 / t oo \ 1 - 

\k R (s,T)\\u(T, R)\ 2 drds 



\0 



< (CR) 1 / 2 yj \u'(s,R)\ 2 ds 
Using Fubini's theorem, the causality condition (ii) and condition (iii), we have 

t oo oo t t t t 

J J \k R (s,T)\\u(T, R)\ 2 drds = J J \k R (s,T)\ds\u(T, R)\ 2 dr = J J \k R (s, r)|ds|it(r, R)\ 2 dr<CR J \u(t, R)\ 2 dr. 
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Therefore, 



1/2 



1/2 



J J k R {s,T)u(s,R)u(T,R)dTds < Cr[ J \u'(s,R)\ 2 ds\ ( J \u(s, R)\ 2 ds 
\o / \0 



(A6) 



The estimates ([A5IA6P imply that 
£(t)-£(0) < - 



1/2 



< 



c 2 

4R 2 



C 2 
4R 2 



\u'(s,R)\ 2 ds 



\u(s,R)\ 2 ds- 



\u(s,R)\ 2 ds- 



C_ 
2R 



\0 



D + S 
R 



D + 5 
R 



t \ 1 2 ' 

\u(s,R)\ 2 ds 

t R 







Je(a)da + ^J J \F{s,x)\ 2 dxds 

— oo 

t t R 

J £ ( s ) ds + §jJ J \F{s,x)\ 2 dxds. 



-oo 



Finally, using the Sobolev type estimate 

R 

2 



\u{s,R)\' 



/f 2R 
d x [u(s, x) 2 ]dx — 2 / u(s, x)u'(s, x)dx < — — S (s), 



we obtain the estimate (|A4[) after the application of Gronwall's lemma. 

As an application of our result, consider the massless scalar wave equation on a Schwarzschild background of mass 
M . Performing a decomposition of the scalar field <& into spherical harmonics F £m , 



oo t 



=o m =-e 



the wave equation reads 



= 0, 



In outgoing Eddington-Finkelstein coordinates r = t + r — r«, p = r, this equation assumes the form 



(A7) 



dl-d 2 



1(1+1) 
P 2 



2M 
P 



1 



Monopolar outgoing solutions have the form (omitting the indices 00 on 

oo 

U(x)dx 



<l>{t,r) = U{n-t) + 2M 



( 2MV 

(t + 2r-n+x) 2 ' " \~R~J 



O 



where U is a smooth and bounded function which vanishes if its argument is sufficiently negative. Since 

2 



(dt + dr. )</>(t, r) 



_ _2M f U{r«-t + 2ry)dy | Q (™\ 



(1 + yf 



\ R J ' 



the boundary condition 
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where 

2M / <j>(0,(R.+2Ry-t).,R)dy 

G{t) -~lv J (TT^ ' 

t/2R 

is perfectly absorbing up to order 2M/R. Performing the substitutions x = r* and y — (t — s)/2R equation (|A7[) can 
be brought into the form (|A1[) and the boundary condition (|A8|) into the form (|A2[) . where 

[0, s > t. 

This kernel satisfies the assumptions (ii) and (iii) with C = M/R, and the potential in Eq. (|A7|) is smooth, positive 
and bounded by the constant W := [t(£+ 1) + 1]/(2M) 2 ; hence, it satisfies assumption (i). For linear fluctuations 
of a Schwarzschild background, Eq. (|A7[) has to be replaced by the Regge- Wheeler or Zerilli equation (f2Tj) . and our 
boundary condition T>2 involves third derivatives of the fields. Therefore, it is not of the form (|A2j) . However, it can 
be reduced to a first order system for which an energy estimate similar to the one presented above can be found by 
introducing auxiliary variables as in the method used in the appendix of Ref. [l8[ . 
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